Quantum instability and edge entanglement in the quasi-long-range order 
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We investigate the build-up of quasi-long-range order in the XX chain with transverse magnetic 
field at finite size. As the field is varied, the ground state of the system displays multiple level 
crossings producing a sequence of entanglement jumps. Using the partial fidelity and susceptibility, 
we study the transition to the thermodynamic limit and argue that the topological order can be 
described in terms of kink-antikink pairs and marked by edge spin entanglement. 
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Strongly interacting many-body systems display a va- 
riety of zero temperature Quantum Phase Transitions 
(QPT) [lj. Quantum criticality is a property of the 
ground state of the system, ultimately arising from a 
reshuffling of the system's energy spectrum when control 
parameters are varied. In many cases, this is accompa- 
nied by a symmetry breaking and by the appearance of 
an order parameter indicating that the macroscopic or- 
der is reached, characterized by long range correlations. 
There also exist phase transitions with vanishing local 
order parameter. A typical example is the Bcrczinskii- 
Kosterlitz-Thoulcss (BKT) transition, occurring, e.g. in 
a system of planar classical spins at a finite temperature 
The transition is characterized by a 'distortion' of 
topological nature of the spatial spin configuration, giv- 
ing rise to quasi-long-range order in the system, whose 
correlation functions display a power-law decay In 
the present paper we will be dealing with ground states 
of one dimensional systems displaying topological order. 
The topological order in one dimension is peculiar and it 
can be characterized by a sensitivity of the ground state 
by varying the boundary conditions, giving rise to certain 
solitonic edge states (see for instance (ij). 

The prospect of practical applications, such as quan- 
tum information processing, has led to an intense ac- 
tivity aimed at a direct inspection of the properties of 
the ground state of a given system. Our purpose is to 
characterize ground states with a non trivial topological 
order through their entanglement content To illus- 
trate what we believe are generic features of topologically 
ordered states with vanishing local order parameter, we 
focus on a specific system: the spin 1/2 XX chain in 
external magnetic field. According to the general equiv- 
alence between the classical-d + 1 and quantum-d crit- 
icality 0], this model displays a QPT from a polarized 
to a BKT phase with quasilong range order. For such 
a model, we demonstrate that the quasi-long range or- 



der manifests itself in the formation of 'entangled edge 
states' (with spins at the edge of the chain sharing en- 
tanglement quite differently from bulk spins) and that 
it is characterized at finite size by an instability of the 
ground state determined by a sequence of energy level 
crossings as the magnetic field is variedM]. This gives 
rise to sudden jumps of both the fidelity [8| and the pair- 
wise entanglement (which, thus, behave non- analytically 
even at finite size). As the size of the chain increases, 
the number of crossings grows until they become dense 
within a sharply defined critical region. At the same 
time, their effect weakens, and the behavior of an infinite 
chain can be obtained through a smearing of the finite 
size observables, a procedure that we apply to the par- 
tial state fidelity. Again, however, this holds true only 
for bulk spins, while those at the boundaries show a kind 
of rigidity and a reduced sensitivity to the QPT. 



THE GROUND STATE OF THE XX MODEL 

We consider A" spin 1/2 particles on a line, coupled by 
nearest neighbor XX interaction, with Hamiltonian 
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where the exchange constant J has been taken as the 
energy unit. In the thermodynamic limit, the system 
undergoes a first order transition from a fully polarized 
to a critical phase with quasi-long range order [3|, [t| . 

Assuming open boundaries (with (Jn+i = 0), we em- 
ploy the Jordan- Wigner and Fourier transformations to 
introduce the fermion operators 
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that diagonalizc H: H = J2k=i A k<4 d * + NB1 with 
A fc = -IB + 2 cos [(vrfc)/(iV + 1)]. The 2 W cigenenergies 

and eigenkets are = X^fcLi ^fe a i^ + -^"-^ an d IVO = 

nf =1 (4) a *' l n >, with 4 l) = (&l4 d *l&) e I ' 1 )- Thc 

state |f2) is the fermion vacuum: dfc|0) = OVA:. 

The ground state and its energy vary with B, and 
different ground states can be classified in terms of the 
number of level crossings occurring in the system as B 
changes. Specifically, when B > 1, < for any k 
and, thus, the lowest eigenvalue is obtained by taking 
the state with = 1 Vfc. The ground state energy re- 
mains e° = —NB as long as B > cos[7r/(7V + 1)]. For 
these values of B no level crossing occurs; the ground 
state is \ip®) = YiiLi^ll^) ■ The first crossing occurs 
at B = cos[7r/(iV + 1)] = B x . For cos[2tt/(7V + 1)] < 
B < cos[7r/(iV + 1)] all of the Aj, are negative ex- 
cept for Ai. Thus, the ground state energy is obtained 
by subtracting its positive contribution: e g = 6® — Ai. 
The corresponding eigenstate is \tpg) = di\ijjg). Let- 
ting Bk = cos\kTr/(N + 1)] and defining the fc-th region, 
Bk+i < B < Bk, we can iterate the procedure above 
to find the ground state \tpg) = dkdk-i • • ■ di\ipg) = 
Yl^Lk+i djl^l) and its energy as 

e* = _(JV-2*)i?- 2^008 (-^) (3) 
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FIG. 1: (color online) Comparison of the entanglement shared 
by spins at the beginning and at the center of the chain. (a) 
One-tangle for spin at site I as a function of the magnetic 
field for I = 1 (solid line) and / = 9 (dotted), (b) Nearest- 
neighbor concurrence between spins 1 and 2 (solid line), and 
9 and 10 (dotted). All of the plots are for N = 19. (c) Partial 
state fidelities and (d) Partial state fidelity susceptibilities as a 
function of external magnetic field Bk = cos( ) for N — 30 
(blue), N = 50 (red), N = 100 (yellow) and N = 200 (green) 
at I — (N + l)/2. The fidelity approaches 1 as JV increases. 
The parameters in the figures are dimensionless. 



FINITE SIZE EFFECTS 



where k represents the number of crossings, 1 < k < N. 
For B < Bn, no other intersection occurs and the ground 
state is simply given by the fermion vacuum state. The 
energy crossings are plotted vs B in Fig. [2] (a). 

To investigate the structure of thc ground state and its 
entanglement content, we re- write it in thc spin language, 
using the eigenstates of of. For k = 0, the ground state 
is \ip°) = | T)®^, which is separable. After the 1-st cross- 



ing, it becomes \<p\) = Y4L1 S l (nL=i a ™ a 



where = 



yj2/(N + l)sm [(nkl)/(N +1)}. This is 
an entangled state, given by a symmetric superposi- 
tion of all possible kets with one flipped spin. After 
k level crossings, the ground state is given by \<Pg) = 



nt =1 [iti sf (n„=r <<>\ 



\)® N . Explicitly 



where ■ 
sites h,h, - ■■ 



h<h<—<lk 
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■ ,lk) is the state with flipped spins at 
, Ik , while the amplitudes are given by 



Ep(-i)X (1) < (2) -^ 



where the 

sum extends over the permutation group. At each cross- 
ing point, the ground state jumps discontinuously in the 
spin Hilbert space from one symmetric subspace to an- 
other, orthogonal to the previous one. 



We now study ground state correlation, entanglement 
and fidelity at finite size. Due to the symmetry of 
the overall state, the one-spin reduced density matrix 
is purely diagonal, pi = idiag(l + (erf), 1 — (cf )); while 
for two spins at sites (l,m), one has pi m = a + \ TT)(TT 
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TIXTJ. I + b-\ 1TXIT I + e(| TIXIT 



l + l 4TXU I), witha± = \[l±(af)±(a*J + (afa* m )}, b ± = 
j[l± (of) T (o^) ~ (of<)] and e = \{afa^). Here, the 
local magnetization and transverse correlation are given 
by (erf) = l-g u and (ofcrfj = (l-gi,i)(l-g m>m )-gl m , 
where 
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that depends on the field B through the index k. The 
longitudinal correlation function (crfa^) is [Tol | : 



(of O = 



G 1,1+1 Gi,i+2 
Gi+1,1+1 Gl+1,1+2 



G m -i t l-\-2 ■ ■ ■ G m —x, 
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where Gi^ m = 8^ m — gi, m - This determinant becomes 
of the Toeplitz type in the thermodynamic limit. For 
nearest neighbors, we simply get (crfcrf +1 ) = Gij+x- 
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With these density matrices at hand, we discuss the en- 
tanglement encoded in the state. Entanglement between 
a single spin and the rest of the chain can be measured 
by the one-tangle t; = 1 — (erf) 2 This quantity de- 
pends on the site for finite chain, and, as a function of B, 
it displays jumps at each crossing point Bk- Specifically: 
(i) ti equals one at zero field for every site; (ii) the jumps 
near B = ±1 become higher and higher moving from the 
end points towards the center of the chain, see Fig. [T] 
(a). Thus, at the onset of the critical region, bulk spins 
are more entangled than end ones. 

To evaluate pairwisc entanglement, we use the concur- 
jo, |e| - y/a+a-Y 
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Its behav- 



Ci >m = 2 max |0, |e| - ^/a 
ior as a function of B for nearest neighboring spins with 
jumps at crossing points is shown in Fig. [T](b). In partic- 
ular, around B = 0, the pairwise entanglement is bigger 
near the end points of the chain; while the reverse occurs 
at the border of the critical region. 

The QPT can be further analyzed through the quan- 
tum fidelity of ground states with slightly different fields, 
B and B = B + SB. Specifically, we consider the partial 
state fidelity of reduced density matrix p a (B) = Tr&p(i?) 
when the system is partitioned as a + b: F a (B 1 B) = 

Tta/ T S /~p a {B)p a {B)\/ p a {B~). F characterizes the degree 
of change of the state as the field is varied. For our 
purposes, it is sufficient to consider the subsystem a to 
consist of just the i-th spinl2j. Within the critical re- 
gion, Fi is unit everywhere except for a series of discrete 
and sharp drops at the crossing points B = Bk- For 
large system sizes, when the crossings are dense in the 
interval \B\ < 1, we can perform a coarse- graining and 
evaluate F only at these points. As a result of this pro- 
cedure, only the drop at B = 1 remains, while all of the 
intermediate ones are smeared out, see Fig. [TJc). Inter- 
estingly enough, this behavior occurs only for bulk spins. 
If, instead, one of the end spins is singled out, the coarse 
grained fidelity stays flat. 

An even more direct evidence that for large N the 
state of the bulk spins changes essentially in a contin- 
uous way except for the critical point B = 1, is ob- 
tained by looking at the fidelity susceptibility [l2|: Xf = 
limsB^o -2\nFi/(5B) 2 , plotted for a bulk spin in FigQ] 
(b). This behavior appears related to the essential sin- 
gularity shown by the block entropy at B = 1 [l3| ■ 



INFINITE SPIN CHAIN 

For TV — » oo, the intervals Bk+i < B < Bk become 
infinitcsimally small and u = (irk)/(N+ 1) becomes con- 
tinuous, lu « arccos(B), so that the sum in ^ gives an 
energy per spin 
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FIG. 2: (color online) (a) The ground state energy per spin 
(black curve) against the magnetic field B in the thermo- 
dynamic limit. The energy for N = 20 is also plotted, to- 
gether with the intersections of the eigenvalues occurring at 
the crossing points Bk = cos[fc7r/(iV + 1)] (vertical lines), (b) 
Thermodynamic limit concurrence between two spins at dis- 
tance r = 1 (blue), r = 2 (red), r = 3 (yellow), r = 4 (green) 
and r = 5 (cobalt). Entanglement between any two spins 
decreases as the distance between the spins increases. It is 
notable that all the parameters in the figures are dimension- 
less. 



which is analytic everywhere within the critical region, 
except for B = ±1. From the finite size analysis, how- 
ever, we know that such region consists of dense set of 
crossings points (see Fig. [5] (a)) and therefore can be 
considered as a line of continuous QPT, with the ground 
state driven by B through various symmetric spin sub- 
spaces with k « (N + l)(arccosi?)/7r flipped spins. The 
correlation functions behave differently for spins close to 
the chain boundaries and for bulk spins. Setting the dis- 
tance r = m — I, the limit N — > oo leads to different 
behaviors of {cr\^a^ +r ) depending on whether I 3> r (bulk 
spins) ori<r (end point spins), [Til ]. In particular, for 
B = (w = tt/2), (af af+J = (l/r 2 )(Z + r) 2 /(l + r/2) 2 , 
with a long distance behavior independent of boundary 
conditions. In contrast, longitudinal correlations are sen- 
sitive to the boundary 



'l u l+r, 



V^A 2 ^ 1 / 2 I » r 

4#(^±i)(V2)i 2 r- 3 / 4 I < r 
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where A = 0.6450025 and K{x) is a function given in[14 
For ./V — > oo, we find gi,i+ r — 2/7r[sinc<j?'/r — sinu;(2^ ■ 
r)/(2l + r)], so that, for bulk spins at B = 0, (erf cj 
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r-2j 
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where [x] is the closest integer 



and (afaf +r ) = [1 - (-1)1 



(-2) 



which 



agrees with 15 1 . 

The concurrence Cij+ r can be derived from the above 
correlation functions. We find that lin^oo Cij +r dis- 
appears for r > 2. Only two nearest spins are entan- 
gled, with liim^ooCu+i = 0.339 at B = 016]. Bulk 
concurrence lim/^oo Ci y i +r between two spins at distance 
1 < r < 5 is plotted for various values of B in Fig. [2] 
(b) , where the decay of entanglement with the distance 
is also shown. For nearest neighboring spins, the con- 
currence goes ~ 1/5, disappearing at B = 1. Near this 
point, bipartite entanglement appears for every r and any 
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two spins in the chain become entangled, although the 
magnitudes become smaller and smaller with r. This is 
because B = 1 is the factorizing point for the XX model, 
with diverging entanglement range 17J. The one-tangle 
limz-joo 77 behaves in a similar way. 



DISCUSSIONS 

In this section we show how the ground state insta- 
bility and the edge entanglement are strictly related to 
the emergence of quasi long range order. It is conve- 
nient to employ a dual basis to describe the system: 
Mn = rim<n °m ■ Once applied to a fully polarized state, 
fi n creates a topological excitation (a kink) [(I, [l^]. In- 
deed, any state with k spin flips (i.e., after k crossings) 
can be viewed as suitable combination of k kink-antikink 
pairs. When B > 1, there are no kinks and the state is 
separable. Near the critical point, with B = 1 — e, the 
ground state consists of a superposition of states with 
a single spin-flip, or a sea of condensed kink-antikink 
pairs (of infinite length for N — * 00); such a conden- 
sation gives rise to the divergence of the concurrence 
range. By decreasing B, the size of the kink-antikink 
pairs decreases and their number increases. At B = 0, 
the ground state has a single (degenerate) kink, with half 
of the spins pointing down and half pointing up. The 
state is 'highly' symmetric and every spin is maximally 
entangled with the rest of the chain, but bipartite entan- 
glement is present with the nearest neighbor only. This 
is due to the fact that the concurrence depends (and it 
is always smaller than) the longitudinal correlation func- 
tion, which ultimately tends to zero because of the pres- 
ence of the kinks. Indeed, with quasi-long range order 
(which arises because the spinwaves are massless), the 
long range correlation function decays since the kinks are 
heavy. The critical region is an instability line because 
the system is driven through different Hilbcrt space sec- 
tors labelled by different quantum numbers (the eigen- 
values of the total magnetization with different 
number of kinks- antikinks pairs. At finite size the phe- 
nomenon of switching among states with a different num- 
ber of kinks is witnessed by the sequence of jumps in the 
entanglement, which are smeared out for N — > 00, except 
at \B\ = 1 as evidenced by the partial state fidelity. 

Kinks are essentially bulk excitations and the picture 
above is modified near the boundaries. Indeed, surface 
spins share entanglement differently from bulk ones: for 
small magnetic fields, every spin is highly entangled with 
the rest, but the end-spins participate essentially to bi- 
partite entanglement, while the bulk ones are rather in- 
volved in multi-partite correlations. On the other hand, 
near \B\ = 1, surface spins are less entangled than those 
at the center (indeed, the first jumps of t\ and Ci 2 
are smaller than the corresponding ones for bulk spins). 
Thus, end spins are less sensitive to the onset of the crit- 



ical region (as shown also by the fidelity), while they are 
more entangled when the quasi long range order is fully 
established. This peculiar edge entanglement is reminis- 
cent of the solitonic edge states found in Q and could 
constitute a fingerprint for topological order in one di- 
mension, valid beyond our specific model. Finally we 
note that preliminary analysis have shown that the for- 
malism we have developed opens the way for the under- 
standing the role of cluster type states [19( close to QPT, 
especially for quantum information purposes. Indeed, a 
Hamiltonian for the cluster state can be mapped onto 
a model with Ising order by means of the dual trans- 
formations employed above. This provides a non-trivial 
disentangling protocol for the cluster states that will 
be explored in a subsequent work. 

Summarizing, through the study of the entanglement 
content and of the fidelity of the ground state, we have 
investigated the quasi long range order in the XX model 
and its connection to the quantum instability arising at 
finite size as a result of the presence of topological exci- 
tations. We have also discussed the special entanglement 
properties of the edge spins arguing that they are a direct 
manifestation of the topological character of the QPT. 
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